THE BENJAMIN-ONO HIERARCHY WITH ASYMPTOTICALLY REFLECTIONLESS 
INITIAL DATA IN THE ZERO-DISPERSION LIMIT 
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Abstract. We study the Benjamin-Ono hierarchy with positive initial data of a general type, in the limit 
when the dispersion parameter tends to zero. We establish simple formulae for the limits (in appropriate weak 
or distributional senses) of an infinite family of simultaneously conserved densities in terms of alternating 
sums of branches of solutions of the inviscid Burgers hierarchy. 



1. Introduction 



The Benjamin-Ono equation is a nonlinear evolution equation governing certain types of internal waves. 
Internal waves are disturbances — set into motion by gravity of the interface between two immiscible 
fluids of different densities. A number of assumptions are employed to deduce the Benjamin-Ono equation 
as a simplified model from the full equations of three-dimensional fiuid mechanics: 

• One supposes that both fluid layers consist of inviscid and incompressible fluids, and (for stability) 
the less-dense fluid rests on top of the denser fluid. 

• One supposes that the waves are propagating in one direction only, which reduces the problem to 
that of two-dimensional fluid mechanics (only the vertical direction and the propagation direction 

survive). 

• One supposes that the top layer containing the less-dense fluid is thin compared to a typical wave- 
length of the interface. This allows the top layer to be treated by a depth-averaging approach. The 
top of this layer is idealized to a rigid horizontal lid. 

• One supposes that the bottom layer containing the denser fluid is infinitely thick. This simplifies 
the contribution from this layer to the dispersion relation of small-amplitude waves on the interface 
in the linear approximation. 

• One supposes that the amplitude of the waves is small compared to the thickness of the top layer. 
This allows the nonlinear effects to be brought in perturbatively, and also means we are assuming 
the interface does not breach the surface by meeting the rigid lid. 

In the linear and dispersionless approximation, deformations of the interface satisfy the one-dimensional 
wave equation. If one perturbatively introduces the balanced effects of weak nonlinearity and dispersion on 
a solution of the wave equation propagating to the right (say) at constant velocity, then the Benjamin-Ono 
equation arises as the first correction in the moving frame of reference of the wave, describing the slow 
variation of the dimensionless wave height u as a function of a spatial coordinate x in the propagation 
direction and the time t: 



Here, subscripts denote partial derivatives, and H is the Hilbert transform operator defined on Z/^(]R) by the 
singular integral 



and e > is a dimensionless measure of the relative strength of dispersive effects compared with nonlinear 
effects. If an initial condition u = uo{x) is given at t = where uq is independent of e, then of course the 



(1) 



ut + 2uux + eUluxx] = 0. 



(2) 
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corresponding solution of ([T]) will depend on e, but we will often not be explicit about this dependence in 
our notation. 

It is quite useful in applications to have accurate and easily analyzable models for internal waves. Indeed, 
one application that is particularly timely is the modeling of submerged "plumes" of oil as were reported 
following the Deepwater Horizon leak in the Gulf of Mexico in May-August 2010. Recent experiments fT] 
performed by Roberto Camassa and Richard McLaughlin at the University of North Carolina, Chapel Hill 
have demonstrated that if oil is emitted as a turbulent jet from an ocean floor leak and a density stratification 
is present in the surrounding fluid, then most of the oil will become trapped at the interface between the 
dense and less-dense fluid layers rather than floating to the surface, even though the oil is less dense still 
than the upper layer. In these circumstances, modeling the motion of the submerged oil plumes within the 
fluid column amounts to modeling the motion of the density interface, that is, modeling internal waves. 

The purpose of this short paper is to show how methods we have recently developed [3] to study the 
asymptotic behavior of solutions of the Benjamin-Ono equation as e J, in the weak topology extend both to 
the whole hierarchy of "higher-order" Benjamin-Ono equations and also to the whole hierarchy of densities 
coming from conservation laws. It is a pleasure to be able to contribute to this special volume of papers 
in honor of Dave Levermore. It will be clear to all of those who have followed his work that our small 
contribution is directly inspired by his groundbreaking analysis with Peter Lax of the zero-dispersion limit 
for the Korteweg-de Vries equation [7j, a project that has had a tremendous impact on fields of study ranging 
from asymptotic analysis of nonlinear evolution equations to the theory of orthogonal polynomials and of 
random matrices. 



2. The Benjamin-Ono Hierarchy 
Let us take the phase space of fields u to be 

(3) ^ ■.= {ue L'^iM.) n C°°(M), uC^) e L^{R), Vfc}. 

This is clearly a linear space over M, and it is also an algebra that is closed under differentiation and Hilbert 
transforms. We will have use below for the Cauchy operators C± densely defined on L^(E) by singular 
integrals as follows: 

(4) C±[/](.x):=lim-^ / JM^ 

The Cauchy operators are bounded with respect to the (M) norm and hence extend to bounded operators 
with the same norms on all of i^(IR). In fact, the operators C+ and — C_ are just the self-adjoint orthogonal 
projections from L^(IR) onto the Hardy subspaces of functions analytic in the upper and lower half-planes, 
respectively. They satisfy the identities 

C+ o (-C_) = (-C_) oC+ = 0, Cl=C+, (-C_)2 = -C_, C+-C_ = l, C++C-^-in. 

Since C± are self-adjoint, it follows that H is skew-adjoint, that is, 

(5) jjn[g]dx^- jjn[f]dx 

whenever / and g area real- valued functions in the phase space . The operators C± and % commute 
with differentiation in x. Also, for functionals I[u\ defined on the phase space we define the variational 
derivative SI/6u by 



(6) j^I[u + tv] 



SI 

— [w] {x)v{x) dx. 
ou 



For functionals / with /[O] = 0, assuming the existence of the variational derivative of / for each u € ^ we 
can recover the functional from its derivative by the formula 

(7) = I f ^[tu]{x)u{x)dxdt^ I (j y[tu]{x)dt\u{x)dx. 



2.1. Conservation Laws for the Benjamin-Ono Equation. The Benjamin-Ono equation conserves an 
infinite number of functionals of u. Tliese may be obtained by several different methods, several of which 
we review (in historical order of discovery) for the reader's convenience. 



2.1.1. The Nakamura Scheme. A. Nakamura [TU] (see also [S]) was the first to deduce an infinite number of 
conserved quantities for the Benjamin-Ono equation ([T]). His derivation is based on a Backlund transforma- 
tion for ([T]). The Backlund transformation is 

(8) — ie/iC+[na;] + 1 — e^" = 

where is an arbitrary parameter. From this equation it can be shown that regardless of the value of ^, 



(9) 



d_ 
dt 



when u satisfies ([T]). Therefore, by expanding n in a power series in jj, with coefhcients depending on u, the 
coefficients will all be densities of conserved functionals of u. Writing n = /irii + /i^n2 -I- /i'^na one easily 
obtains a recurrence in which the conserved density n„ is explicitly given in terms of ni, n2, . . . , rin-i and 
u. In fact, the recurrence can be made more explicit by noting that the nonlinearity of the scheme is only 
quadratic in Ux'- differentiating ([s]) with respect to x and then using ([8| to eliminate e^" from the result 
one arrives at 



(10) 

and therefore ni ; 



Tlx = fJ. (uux + Ux + ienxC+[nx\ + ieC+[nxx]) , 



and 



(11) 



m-2 

Unm^l,x + i<^C+[nra-l,xx] + X! "-i.aiC+Kn-l-j.a;], m = 2,3,4,.. 



The first several densities are given by: 



(12) 



ni = u 
1 
2 
1 



n2 



"3 
714 



u'^ + e 



+ e 



-uHlux] + iuti 



1 



-T-Lluux 



-i'HWx 



1 



-uniuux] 



2I 



e -iu^ + + eM -MH[ux\ - -uux + -H[uH[ux\_ 



e -UxHlUx 



4 ' ^' 4 



+ e'^^i'H [{uxf + uuxx] ■ 

In the expression for 71,4, only the terms on the first line contribute to the integral over K. Indeed, those on 
the second line are derivatives of functions in and those on the third line have zero integral because H is 
skew-adjoint and = — 1. To deduce that the terms on the fourth line have zero integral, it is necessary to 
correctly interpret the integral as the integrand is not of class L^(M). The key identity here is the following: 
if / e ^ then 



(13) 



lim 

-Rtoo J — R 



H[f]{x) dx = 0. 



2.1.2. The Fokas-Fuchssteiner Scheme. The following method is due to Fokas and Fuchssteiner [5]. It is 
based on Lie-theoretic analysis of one-parameter symmetry groups of the Benjamin-Ono equation Q. One 
begins with 



(14) 

and then recursively defines 
(15) 



6Im 

5u 



1 <5 



.-2 5u 



m = 1 and ^—\u\=u 
ou ou 



2xuUx -f u + e xTi, \uxx\ H — T~L [ux 



SIm-l 

Su 



[u] dx, m = 3, 4, 5, . . 



To compute the variational derivative on the right-hand side one needs the identity 



(16) 



H[x,f]=xH[f] if / f{x)dx^Q. 



Although the explicit function x appears in the integrand on the right-hand side of (15), the recursion 



guarantees that the variational derivatives produced are all generated from u and its derivatives and a finite 
number of applications of %. For example, 



(17) 



5h 

Su 



S_ 
Su 



2xu^Ux + u^ + e \ xuT-L[uxx\ + 7:uH[ux\ 



dx 



el-Llux] 



and 



(18) 



Sli , . Q / 3 . 



-Hluux 



As X does not appear explicitly in the resulting expressions, it follows from the formula Q that the functionals 
Ijn are integrals of densities that also do not involve x explicitly. These densities may be easily obtained from 
the variational derivatives simply by first multiplying through term-by-term by u and then dividing each 
term by its homogeneous degree in u. For example, from each SIm/Su we obtain a corresponding conserved 
density /,„ as follows: 



(19) 



/i = u 

h = + £^un[ux] 

h = ^w"* + e (^u'^'H[ux] + ^un[uux 



2I 

e -uUr, 



Unlike the densities produced by the Nakamura scheme, these densities are all absolutely integrable if u € 



2.1.3. The Kaup-Matsuno Scheme. The following method was derived from the inverse-scattering transform 
for the Benjamin-Ono equation by Kaup and Matsuno [5]. Set ki := u and then define recursively 

(20) km := uC+[kjn-i] + ie ^ ^"^ , m = 2,3, 
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Then, the quantities km are all densities of functionals conserved by ([T]) . The first several densities obtained 
by the Kaup-Matsuno scheme are 
ki = u 

k2 = -u^ — -iuHlu] 

ks — -u^ — -u7i[u7i[u]] — -iu^Hlu] — -iuHlu^] + e ^-iuu^ + -uH[ux 

fc4 = lu^ - lu^niuniu]] - luniu^niu]] - lun[un[u^]] + luniuniuniu]]] 

(^J-j 8 8 8 8 8 

8 8 8 

/I 3 1 3 1 1 

+ e I -u^T-LluJ + -uHlnuJ + -uUxHlu] + -iuux — -iuHluT-LluJ] — -iuT-LluxT-Llu]] 
\2 44 4 2 4 

As is the case for the Fokas-Fuchssteiner scheme, and unlike the Nakamura scheme, the densities produced 
by the Kaup-Matsuno scheme are all absolutely integrable for m £ 3^. 

2.1.4. Comparison of the Schemes. Although the three schemes clearly do not give rise to identical densities, 
they apparently produce exactly the same integrals. That is, for w e <f3^ we typically have rim, fm, and km 
being different (unequal) expressions generated by derivatives of u and applications of H. However, 

(22) Im[u] -.^ / Umdx^ / fmdx^ kmdx, m= 1,2,3, 

Jr Jr Jr 

Here in the case of rim dx the integral must generally be interpreted in the "principal value at infinity" 
sense as in (13). From each of the three schemes it follows that if m is a smooth solution of ([I]), then 



dim/dt — for all m. We do not provide a direct proof of the equivalence of the integrals generated by each 
of the three schemes here, although in the case where e = all three schemes produce the same result: if 
u G ^ is independent of e, then 

(23) Mmlmit) ^ — [ u(x)"'dx. 

This is rather straightforward to show from the Nakamura and Fokas-Fuchssteiner schemes, while for the 
Kaup-Matsuno scheme it follows from a lemma proved in the appendix of [5]. In any case, this fact easily 
establishes the functional independence of the integrals Im[u], m = 1, 2, 3, ... . 

The variational derivatives SIm/Su in fact play a dual role, as is proved in both [5] and [8j (in the respective 
context of two different schemes): they also serve as densities of integrals: 

(24) [ — ^ dx — {m~ l)Im-i[u\, TO > 2, or Im[u] ~ [ dx, m > I. 

Jr ou to du 

The main advantage of the Nakamura scheme is that it can be used to place all of the equations of the 
Benjamin-Ono hierarchy (see below) in bilinear form after which Hirota's method can be applied to deduce 
the form of the simultaneous A^-soliton solution of the entire hierarchy. This is done in [5] , and we will use 
the resulting formulae below. However, due to the presence of terms of the form T-L [/] for / £ that are 
not themselves of the form f = ioi g G ^ , the densities generated by the Nakamura scheme are generally 
not absolutely integrable and require a more careful interpretation of the integral. 

The advantages of the Fokas-Fuchssteiner scheme are (i) that it provides a direct one-term recurrence for 
the variational derivatives of the conserved quantities (which unlike densities are uniquely determined by 
the functional), (ii) that the variational derivatives it generates are manifestly real and the corresponding 
densities fm are simpler than in either of the other two schemes, and (iii) that it allows a direct proof of the 
fact that the functionals Im are all in involution with respect to the Poisson bracket: 

(25) {/,J}:= / ^(x)^^(x)dx, 

./b du ox ou 



in other words, {Ij,Ik} = for all j,k. 

In the Kaup-Matsuno scheme, the quantities Nm '■= km/u are the coefficients in the Laurent expansion 
about A = oo of the eigenfunction N{X;x,t) that is a simultaneous solution of the two linear equations 
making up the Lax pair for the Benjamin-Ono equation. Here A G C is the spectral parameter. Because the 
integrals I,n[u] :— Jj^kndx are obtained by expanding a scattering eigenfunction, they equivalently encode 
the scattering data in an explicit way, and in 6J one can find explicit formulae for /,„ in terms of the 
discrete spectrum and the reflection coefficient of the direct scattering problem. Thus, the Kaup-Matsuno 
scheme leads not just to an infinite collection of integrals of motion of the Benjamin-Ono equation, but also 
a hierarchy of trace formulae equivalently expressing each both as a functional of u and also as functional 
of the scattering data. These formulae can be used to deduce from initial data the asymptotic distribution 
of eigenvalues Xj of the scattering problem in the zero-dispersion limit (see and section 3.2 of [9j). 

2.2. Construction of the Hierarchy. The Benjamin-Ono equation ([T]) can be written in Hamiltonian 
form as 

(26) ut = - — — . 

ox ou 

The Noetherian symmetry of ([!]) associated with the conserved quantity /„,, is the Hamiltonian flow with 
Hamiltonian /„: 

(27) ^,^=_^^, fc = l,2,3,.... 

ax ou 

Here tk is the parameter of the symmetry group generated by Ik+2- The fact that the integrals are all in 
involution 5J implies that these flows are all compatible (that is, the symmetry group is abelian), so given 
a smooth function uq and a positive integer K, there will exist a function u{x, ti,t2, ■ ■ ■ ,tK) satisfying 



u{x^ 0, 0, . . . , 0) — Uo{x) and equations (27l for A; = 1,2, . . . , K . Equations (27) constitute the Benjamin-Ono 
hierarchy. 

3. Setting up the Zero-Dispersion Limit 

3.1. Formulation of the Problem. The problem we wish to consider is the following. Let uq G ^ he 
given, an initial condition independent of e. For each e > we may construct the simultaneous solution 



u{x,ti,t2, ■ ■ ■ ,tK) of the Benjamin-Ono hierarchy (27) of commuting flows satisfying the initial condition 



u(x,0,0, . . . ,0) = uo{x). The question of interest is the asymptotic behavior of u(x,ti,t2, ■ ■ ■ ,tK) in the 
zero-dispersion limit e J, 0. As a ffist step, we will address this problem by establishing the existence of the 



dispersionless limits (in appropriately weak topologies) of all of the conserved densities (see (24)) 

(28) D,,,:^-^-^, m = 1,2,3,.... 

TO Ou 

Note that Di — u. In general, Dm differs from all three of rim, /m, and km by a "trivial" density that 
integrates to zero for all u £ However, the densities Dm are those that most directly yield a dispersionless 
representation. 

3.2. Admissible Initial Conditions. We will further assume that uq G ^ satisfies the following conditions 
adapted from [9 : 

• uo(x) > for all a; e M. 

• There is a unique critical point xq G M for which Uq{xq) — 0, and mJ^xq) < making xq the global, 
nondegenerate maximizer of mq- 

• Uq exhibits power-law decay in its tails: lim2,_j.-|-oo "0(2;) — and 

(29) lim |a;|«+^M^(a;) = C± for some q>\, 

where C+ < and C_ > are constants. 

• For each fc = 1, 2, . . . , let /(x) = u^ixf . The n in each bounded interval there exist at most 
finitely many points a; = ^ at which /"(C) = 0, and each is a simple inflection point: /"'{£,) 7^ 0. 

Such uq G ^ will be called admissible initial conditions. An example is shown in Figure [l] 

6 



Figure 1. An admissible initial condition for which Uq{x) = — 5x(l + cos(7ra::)/2)/(H-a;^)^. 
Left: Uo{x). Center: Uq{x). Right: Uq{x). Note that this particular admissible initial 
condition has an infinite number of inflection points asymptotically near integer values of x 
for X large. 



3.3. Reflectionless Modification of the Initial Data. Our method will be to study the Benjamin-Ono 
hierarchy with admissible initial condition uq using the inverse-scattering transform for the x-part of the 
Lax pair (see IH IHI)- The first step in the process is to associate to uq its scattering data, consisting of a 
complex- valued refiection coefficient /3(A), A > 0, as well as discrete eigenvalues {Aj < 0} and corresponding 
phase constants {7^ € M}. Matsuno was the first to observe that the conservation laws can be used to 
deduce information about the scattering data (see section 3.3) when the parameter e > (which appears 
parametrically in the scattering problem although the potential uq is independent of e) is small. His analysis 
becomes more rigorous with the use of the trace formulae arising from the Kaup-Matsuno scheme. For 
admissible initial conditions, Matsuno's main results are: 

• The reflection coefficient /3(A) is small when e <C 1. 

• The number iV of eigenvalues is large when e <C 1 , but 



(30) WmeN = M := / uo{x)dx. 

The number N[a, b] of eigenvalues in the interval [a, 6], — L < a < b < 0, L :— maxxes. uai^), satisfies 

(31) \imeN[a,b]= f F{X)dX, F(A) := -!-(a;+(A) - a:_(A)). 

Here a;_(A) < a;+(A) are defined for —L < A < as the two solutions of the equation uo{x) = — A. 
They play the role of turning points in this theory. 



(Note that the "mass" M defined by ( 30 1 is finite for admissible uq although it is not so for general elements of 
The solution of the Benjamin-Ono hierarchy for an admissible initial condition is therefore in particular 
approximately reflectionless in the zero-dispersion limit. In the absence of reflection the exact solution of 
the hierarchy is a multi-soliton solution that takes the form [SI: 

(32) u{x,ti,t2, ...,tK)^ 2e — 3{log(r(a;,ii,t2, • ■ ■ ,tk))} 
where the "tau function" is 

(33) t(x, ^1, ^2, ■ • ■ , tK) ■■= det(I + ie-^A) 
and A is an X iV Hermitian matrix with constant off-diagonal elements 



(34) Anm ■= ^ n^m 
and diagonal elements depending explicitly on x,ti,t2, ■ ■ ■ , tx'- 

(35) Ann ■= -2A„ {X{Xn;X,ti,t2, . . . ,tK) + Jn) , 

where for future convenience we define a polynomial in A by 

K 

(36) X{X; x,h,t2,..., tK) -.^x- ^(fc + l){~X)Hi,. 

k=l 



To specify an appropriate family of exact solutions of the Bcnjamin-Ono hierarchy, first define the exact 
number of approximate eigenvalues by 

M 



(37) 



N{e) 



Then, define approximations {~\n\'^=i with — L < Ai < A2 < • ■ ■ < Ajv(£) < by quantizing the Matsuno 
eigenvalue density: 



(38) 



y" V(A)dA = e(^n-0, n = f,2, 



,7V(6). 



7(A„), ^{X):=--{x+{\)+x-{\)), -L<\<0. 



Next, define approximations for the phase constants {7,1} by setting 
(39) 7„ ' ^ 

Now, for each e > 0, let u = u{x,ti,t2, ■ ■ ■ ,tK) denote the exact solution of the Benjamin-Ono hierarchy 
given by the reflectionless solution formula (32 ) with determinantal tau function f involving the N{e) x N{e) 
Hermitian matrix A whose elements are given by ( 34 1-( 35 ) with A„ and 7„ replaced by A„ and 7„ respectively, 
for 1 < n < N{e). In [3] it is proved that for admissible uq, 



(40) 



lini 



\uo{x) - u{x, 0, 0, . . . , 0)p = 0, 



so that the replacement of the scattering data we have just made amounts to a modification of the initial 
condition that is negligible in the L^(K) sense in the zero-dispersion limit. 



4. Distributional Limits of Conserved Densities 

The conserved densities Dm , m > 1 (we are using tildes to remind the reader that these are expressions in 
u, its derivatives in and Hilbert transforms thereof), all have representations in terms of the tau function 
associated with u. Indeed, since u satisfies ([27|) it follows from using the formula (32) for u that 

(41) 



Integration in x using decay at a; = ±cx) to fix the integration constant yields the formulae 

2e d 



(42) 



Dr. 



-5{log(f)}, m = 2,3,...,X + l. 



Of course since Di = u a. slightly different formula holds for Di according to (32). In principle, this gives 
a way of evaluating Dm for arbitrary m, although one must include dependence on a sufficient number of 
times tk by choosing K large enough. 

Let ai < a2 < ■ ■ ■ < ctN{e) denote the (real) eigenvalues of A. Then, we may write the densities in the 
form 



(43) 
where 

(44) 



u = Di = 



dU 
dx ' 



Dm = 



1 dU 

m dtm-i 



m = 2,3,...,i^ + l. 



U{x,ti,t2, ■ ■ ■ , Ik) '■= e 2 arctan(e ^a^). 



Remarkably, the function U has a completely explicit zero-dispersion limit: 
Proposition 1. Uniformly on compact subsets ofM.^^^, 



(45) 



lim U(x, ti,t2, 



,iK) = V{x,ti,t2 



,tK) 



TT sgn(a)G(Q;; x, ti , t2 



Ik) da 



wh 



ere 



1 r° . .dx 



(46) G{a;x,ti,t2,...,tK) := -— J X/(a) ^ 
and where X/(q^) denotes the indicator function of the interval 

(47) / := [-2A {X{X; xMM,---,tK)~ x+{X)) , -2A {X{X; x^M, ■ ■ ■ ,tK) - (A))] . 

Proof. This is a simple generalization ol Proposition 4.2 Irom 9J and it is proved in exactly the same way 
(see in particular sections 4.1-4.3 ol that reference). For the reader's convenience we will simply describe 
the idea ol the prool. 

The key observation is that the eigenvalues ol the matrix A have a limiting density G{a] x,ti,t2, . ■ . , tx)] 
that is, the normalized (to total mass Af) counting measures ol eigenvalues ol A converge in the weak-* 
sense to G da as e | 0. This lact is proved using Wigner's method ol moments. One studies the asymptotic 
behavior ol traces ol arbitrary powers ol the iV(e) x N{e) matrix A in the limit N{e) — >■ oo, and with 
the use ol some combinatorial arguments and approximation in terms ol diagonal and Toeplitz matrices 
one obtains leading-order asymptotic formulae lor these, which in turn are proportional to moments ol the 
eigenvalue counting measures. Then one solves the moment problem lor the limiting moments to obtain 
Gda. Finally, by estimating the extreme eigenvalues of A one is able to convert convergence ol moments to 
weak-* convergence. 



Next, one observes that the exact formula (44 ) can be written as the integral ol the lunction 2 arctan(e^^a) 
against the normalized counting measure ol eigenvalues ol A. Pointwise, the integrand converges to tt sgn(Q;) 
as e 4, 0, and by a carelul dominated convergence argument one then establishes the desired locally uniform 
convergence ol U to V. □ 

Before giving our next result, we recall the inviscid Burgers hierarchy. Consider the equation 
(48) UB = uo(^x- ^(fc + l)u|<fc^ . 

For ti , . . . , all sufficiently small (given a; e M) it follows Irom the implicit lunction theorem that there 
exists a unique solution UB{x,ti,t2, ■ ■ ■ ,tK) ~ uo{x). As tk increases Irom zero, there will be bilurcation 



points at which the number ol solutions ol ( 48 ) increases by a finite even integer (this is due to the condition 
on inflection points ol f{x) = uq{x)'' satisfied by admissible uq). Therefore, near a given x and at given 
values ol ti,t2, ■ ■ ■ ,tK , there will generically be an odd finite number 2P{x, ti,t2, ■ ■ ■ , tx) 4- 1 ol distinct 



solutions mb.o < ^*b.i < • • • < ub,2P to (48), and each is differentiable with respect to all of the independent 



variables x and ti,t2, ■ ■ ■ ,tK- By differentiation of (48) one observes that each of the solution branches is a 



function UB{x,ti,t2, ■ ■ ■ jtx) that simultaneously satisfies the equations 



These are the partial differential equations of the inviscid Burgers hierarchy. The simultaneous solution of 
these equations with initial condition ub{x, 0, 0, . . . , 0) = ^0(2:) is accomplished by the method of character- 



istics and produces the solution in implicit form (48 1. We note that whereas typically when Burgers-type 



equations appear in the theory of partial differential equations one is interested in single-valued weak solu- 
tions representing shock waves, our interest here is in the multivalued solution produced by finding all real 
solutions of the implicit equation ( 48 1 . See Figures [l]and[2) 

Proposition 2. The function V is of class C^{M.^^^), and 

dV 

(50) — = ^(-l)"uB,n(a:, tl, ^2, . . . , tK) 

n=0 

while 

(51) ^-j- ^y^{~ir-UBAx,h,t2,...,tKr, m = 2,3,...,i^ + l, 

mdtrn-1 ^ m 

n—O 



Figure 2. Left and center: the solution ub of (48) with initial data as specified in Figurejl] 
with tk — for fc > 3 and with ^1=^2=^3 = 0.1 and ti = t2 = ^3 = 3 respectively. Right: 
a close-up of the case when <i = <2 = ^3 = 3 displaying intervals of x with P = 0, 1, 2 (one, 
three, and five branches, respectively). 

where P = P(x,ti,t2, . . . ,tfc)- Both of these formulae assume that {x,ti,t2, ■ ■ ■ ,ti<-) G M^'+^ is a point at 
which the integer P is well-defined; however, since new solution branches bifurcate in pairs from the same 
value of Ub it is clear that the formulae extend by continuity to all of M^"*"^ . 

Proof. The proof is virtually identical to that of Lemma 4.12 from 9 . The idea is as follows. With the 
use of the explicit formula for G{a] a;, ti, ^2, ■ • ■ , tx) one can exchange the order of integration in the formula 



(45 1 to obtain 





(52) V{x,ti,t2,...,tK) ^ J J{X;x,ti,t2,-..,tK)dX 
where 

f-7rP(A), X{X;x,ti,t2,...,tK) < X-{X) 

(53) J{X;x,ti,t2,...,tK) := < X(A; a;, ii, ia, ■ • ■ , ) + 7(A), a;-(A) < X{X; x,ti,t2, . . . ,tK) < x+^X) 

[nF{X), X{X;x,ti,t2,...,tK) > x+{X). 

The integrand therefore has a different form as a function of A in three different types of subintervals of 
[— L, 0], with boundary points given by the solutions A of the equations 

(54) X{X;x,ti,t2,...,tK)^x±{X). 

Recalling that a;±(A) are two branches of the inverse function of uq in the sense that for —L < A < 0, 
uo{x±{X)) — —A, both of these equations can be combined in the form 

(55) -X = uoiX{X;x,ti,t2,...,tK)), 



which one immediately notices is the same implicit equation ( 48 ) providing the multivalued solution of the 



Burgers hierarchy, under the substitution ub = —A. Differentiation of (52 1 using Leibniz' rule to take into 



account the moving boundaries then yields the desired formulae. □ 
Now we may formulate our main result. 



Theorem 1. Let v e L^{R). Then 



(56) hm y Di{x,ti,t2,...,tK)v{x)dx^ / \y {-iy'^UB^n{x,ti,t2, ■ ■ ■ .tx) ] v{x) dx 

holds uniformly for {ti, . . . , tx) in compact subsets ofM.^. Recalling that Di = u, this means that u converges 
to the alternating sum of branches of the multivalued solution of the Burgers hierarchy with initial condition 
Uq in the weak L'^{M.x) sense. Moreover, the convergence is in the strong L^(IR^) sense if ti, . . . ,tK are all 
sufficiently small that ub is single-valued as a function of x. 

Now let (j) G 2>{M.) be a test function. Then for m = 2, 3, . . . , X + 1, 



(57) lim/ Djnix,ti,t2,...,tK)4'{tm-l)dtm-i^ [\y^{-l)"—UB,n{x,ti,t2,...,tK)"'\4'{tm-l)dtr. 

'^io Jk Jr \^;^ m I 
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holds uniformly for {x,ti, . . . , t„i, . . . , Ik) in compact subsets of M.^ . Therefore £),„ converges to the 
alternating sum of m-th powers, weighted by 1/m, of the branches of the multivalued solution of the Burgers 
hierarchy with initial condition uq in the topology of ^'(]Rt^_J, that is, distributional convergence with 
respect to t„i-i • 

Proof. The distributional convergence of Dm for m > 2 clearly follows from our above results. It is also 



easy to conclude that (561 holds if v is specialized to a test function e To strengthen this to weak 

L^(IRa;) convergence and strong L'^{Rx) convergence pre-breaking, one follows nearly verbatim the arguments 
on pages 254-256 of ij. □ 

We expect that with some additional effort, the nature of the convergence of D„i for m > 2 can be 
strengthened to exactly the same type as is available for Di = u, a type of convergence that is more suitable 
for evaluation at a point in the phase space 3^ of fields. This expectation is based on the reasonable 
hypothesis that the weak (or distributional) nature of the convergence stems from the presence of wild 
oscillations that can be modeled by modulated P-phase wave exact solutions of the Benjamin-Ono hierarchy 
as have been described by Matsuno [5j using the bilinear method of Hirota. These P-phase waves have also 
been obtained directly from the Lax pair for ([T]) (the k — 1 case of the hierarchy only) by Dobrokhotov 
and Krichever [3], who further provided a formal Whitham-type modulation theory for these waves, noting 
that the modulation equations simply take the form of 2P + 1 copies of the inviscid Burgers equation 



(equation (49 1 for k = 1). In hght of our results, it appears that these 2P + 1 copies should be globally 
viewed as sheets of the same multivalued solution. In any case, if one interprets the distributional limits 
in ^'(Mi^^ j) as local averages of Dm over vanishingly small intervals of tm-i, then assuming only that the 
wavenumbers and frequencies are not rationally dependent, these averages could just as well be calculated 
over small intervals in x, holding ti,t2, ■ ■ ■ ,tK fixed. In other words, if the weak limits are necessary due to 
the presence of modulated multiphase waves of wavelengths and periods proportional to e, then there should 
at generic points be no difference between convergence in &'{M.t^_-^) and convergence in ^'(M^). A proof 
of such a result probably requires resolution of the microstructure as could be obtainable from an approach 
to the zero-dispersion limit that starts with the nonlocal Riemann-Hilbert problem of inverse scattering 
for Benjamin-Ono, and that involves the development of some new analogue of the Deift-Zhou asymptotic 
method as has been applied [2 to strengthen the zero-dispersion limit of Korteweg-de Vries equation. We 
hope to be able to announce progress in this direction in the near future. 

It seems to us that while the Benjamin-Ono equation looks at first glance to be a more complicated model 
for wave propagation than the more famous Korteweg-de Vries equation due to the presence of the Hilbert 
transform and its concomitant nonlocality (and perhaps even at "second glance" , since the treatment of the 
Benjamin-Ono equation by the inverse-scattering transform method is far less well-understood than in the 
case of the Korteweg-de Vries equation), in fact it is far simpler in the zero-dispersion limit. Indeed, the 
asymptotic formulae that are the analogues in the Korteweg-de Vries case of our limiting formulae for Dm 
require the solution of a variational problem for a quadratic functional with constraints as was found by 
Dave Levermore and Peter Lax in their pioneering work [7 , while for Benjamin-Ono it suffices to be able to 



solve the implicit algebraic equation ( 48 ) for ub , or alternatively to solve the system of partial differential 
equations (49 1 numerically by the method of characteristics. These are far more elementary tasks. We 
want to stress this point to hopefully encourage the use in the practical modeling of internal waves of the 
simple approximate formulae available for the Benjamin-Ono equation and its hierarchy when the dispersion 
parameter e can be reasonably assumed to be small. 
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